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effectiveness of a system prior to installation. This article discusses the five main parametric models for 
photovoltaic (PV) systems available in the literature. The current-voltage characteristics of a photovoltaic 
module can be reproduced modelling the PV panel as an equivalent electrical circuit made of linear and 
Keywords: non-linear components. The parameters describing such components are directly related to the 
Five parameter model performance characteristics of the specific PV panel, which are generally available in a graphic form 
en with respect to standard values of temperature and incident irradiance. Five of the most recent models 
Renewable eiiergy are analysed in detail describing the algorithm, showing how the parameters were determined and 
highlighting the most important mathematical and physical assumptions that characterised each model 
and the calculation process. To simplify the examination and reading of the procedures referred to in the 
five models, all equations are rewritten using the same nomenclature. Furthermore, to assess the 
accuracy of each model with respect to the data provided by the manufacturer, this paper compares the 
I-V and P-V curves at various temperatures and irradiance for a generic PV panel. 
© 2014 Elsevier Ltd. All rights reserved. 
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1. Introduction 


A steady increase in hydropower and the rapid expansion of 
wind and solar power has cemented the position of renewables as 
an indispensable part of the global energy mix; by 2035, renew- 
ables will account for almost one-third of the total electricity 
output. Solar power is growing more rapidly than any other 
renewable technology. Renewables will become the world’s 
second-largest source of power generation by 2015 (roughly half 
that of coal), and, by 2035, they will rival coal as the primary 
source of global electricity |1]. 

In this framework, photovoltaic generation systems have the 
opportunity to be suitable for use due to their important advan- 
tage in being able to produce electrical energy very close to the 
electric load. In this way, transmission losses are avoided, and it is 
also possible to satisfy the peaks in the daily load diagrams 
because they supply the maximum power in correspondence to 
the maximum request. Moreover, photovoltaic plants do not emit 
pollutant emissions or vibrate; due to their modularity, they can 
comply with the morphology of the installation sites [2]. 

Promoting the dissemination of renewable energy sources is 
essential for accurate assessment of power production to estimate 
the economic and environmental performance with a high grade 
of reliability [3,4]. 

In the field of PV systems, it is essential to have a good 
simulation tool that can forecast the physical behaviour of a PV 
panel by varying the solar irradiance and the operating tempera- 
ture. The operating temperature is also a key variable of the 
forecasting problem because it directly affects the electrical out- 
puts (current, voltage and power) and the efficiency of the solar 
device. One of the thermo-physical processes that most affects the 
value of the operating temperature is represented by convective 
and radiative heat exchange between the panel and the surround- 
ing environment. In [5], it was observed that the convective heat 
transfer coefficient is strongly influenced by the wind speed, 
although it is less affected by the wind direction and is almost 
unaffected by the air temperature. Concerning the radiative heat 
exchange, the operating temperature also strongly depends on the 
glass transparency and the emissivity and absorption coefficient of 
the cell. Some authors have attempted to estimate the operating 
regimen of a photovoltaic panel connected to an electrical load in a 
simple and precise model, noting the operating temperature [6-8]. 

Some existing software, such as PVFORM [9,10], PVWATTS [11], 
PHANTASM [12], TRNSYS [13] and PVsyst [14], which is quite 
sophisticated and is intended for advanced users, evaluates the 
physical behaviour of a photovoltaic component and requires a 
large number of parameters, which are difficult to find in technical 
data and which are often published without the correlations 
used. Consequently, normal users usually employ the databases 
included, sometimes without understanding the importance of 
these data in the assessment of the power output. For these 
reasons, it is certainly wise to learn more about the features and 
the mathematical approximations of the various models used to 
describe the thermo-electrical behaviour of PV systems. 

The algorithms and models that can be used to assess the PV 
power output are generally of two types: those that use selected 
points of the I-V curve [15-19] and those that exploit the entire 


characteristic curve [20-23]. The first group of algorithms involves 
the solution of five equations, using thermal drift coefficients, 
measuring the slope of the I-V curve and exploiting selected 
points, i.e., the open-circuit, short-circuit and maximum power 
points. Although calculating the exact solution of these equations 
requires iterative techniques, this method is often much faster and 
simpler compared to the second group of models. Indeed, the 
group of algorithms based on curve fitting, using deterministic 
optimisation algorithms, has a disadvantage that leads to solutions 
for several sets of local minima [24]. 

Within the first group of algorithms, the analytical models that 
are employed the most evaluate a solar cell as an equivalent 
electrical circuit composed of a set of generally non-linear com- 
ponents. The electrical parameters of the equivalent circuit are 
directly related to the performance characteristics of the photo- 
voltaic component, which are generally described in graphical 
form for standard values of temperature and irradiance. 

The paper is organised as follows: Section 2 presents the types 
of PV, related topics and their performances; and Section 3 
presents the equivalent circuit for PV modules and describes the 
five most recent and most cited mathematical models. An outline 
of the parameter extraction procedure is given for each model. 
Section 4 presents a discussion of the results obtained along with 
the procedures used. 


2. Solar photovoltaic technologies 


The photovoltaic technologies are generally divided into five 
categories: 


1. Multi-junction solar cells: Multi-junction solar cells, or tan- 
dem cells, are solar cells containing several p-n junctions. Each 
junction is tuned to a different wavelength of light, which 
reduces one of the largest inherent sources of losses, thereby 
increasing the efficiency. Traditional single-junction cells have 
a maximum theoretical efficiency of 34%; a theoretical “infi- 
nite-junction” cell would improve this efficiency to 87% under 
highly concentrated sunlight. 

2. Single-junction GaAs: Gallium arsenide is an important optical 
compound semiconductor and is a nearly ideal material for 
high efficiency single-junction solar cells. In addition, GaAs is 
capable of withstanding radiation damage for a long time. 
Owing to its high efficiency and high resistance to radiation, 
GaAs cells have been favourites for powering satellites and 
other spacecraft [25]. 

3. Crystalline Si cells: Silicon technology has been the dominant 
one for supplying power modules in photovoltaic applications, 
and the changes likely are an increasing proportion of multi- 
crystalline silicon and mono-crystalline silicon being used for 
high-efficiency solar cells, while the use of thinner wafers and 
ribbon silicon technology continues to grow [26,27]. 

4. Thin-film technologies: Thin-film solar cells are basically thin 
layers of semiconductor materials applied to a solid backing 
material. Thin films greatly reduce the amount of semiconduc- 
tor material required for each cell compared to silicon wafers, 
and hence, lower the cost of production of photovoltaic cells. 
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Nomenclature 

a ideality factor 

D diode diffusion factor (approximately constant) 
G irradiance [W/m?] 

Gsrc irradiance at SRC [W/m?] 

I current generated by the panel [A] 

Ip diode current [A] 

I, photocurrent [A] 

Imp current at the maximum power point [A] 
Isc short circuit current of the panel [A] 
Io reverse saturation current [A] 

K thermal correction factor [Q/°C] 

k Boltzmann constant [1.381 x 107” J/K] 
J thermal correction factor 

M air mass 

m constant 

n diode quality factor 

N; number of cells wired in series 

q electron charge [1.602 x 10719 C] 

Rı electrical load [Q] 

Rs series resistance [Q] 


Rso reciprocal of slope of the I-V characteristic of the 
panel for V=Voc [V] and 1=0 [Q] 

Rsh shunt resistance [Q] 

Rsho reciprocal of the slope of the I-V characteristic of the 
panel for V=Voc [V] and I=0 [Q] 

SRC Standard Reference Conditions 

T cell temperature of the PV panel [K] 

V voltage generated by the PV panel [V] 

Vinp voltage at the maximum power point [V] 

Voc open circuit voltage of the panel [V] 


Greek symbols 


a ratio between the current irradiance and the 
irradiance at SRC 

Eref material band-gap energy [1.12 eV for silicon] 

ESRC amplitude of energy gap [J] 

y shape factor 

Misc thermal coefficient of the short circuit current [A/°C] 

Hv,oc thermal coefficient of the open circuit voltage [V/°C] 


GaAs, copper, cadmium telluride (CdTe), indium diselenide 
(CulnSez) and titanium dioxide (TiO2) are materials that have 
been used the most for thin-film PV cells [27]. 

5. Other emerging photovoltaic: In this category, there are different 
types of PV technology, such as dye-sensitised, organic, inorganic 
and quantum dot solar cells. 


A dye-sensitised solar cell (DSSC, DSC or DYSC) is a low-cost 
solar cell belonging to the group of thin-fOilm solar cells [28]. It is 
based on a semiconductor formed between a photo-sensitised 
anode and an electrolyte that constitutes a photo-electrochemical 
system. 

Organic photovoltaic (OPV or excitonic) cells may be based on 
organic dyes, semiconducting polymers, small semiconductor 
molecules, or on some combination of these species. In all cases, 
however, light absorption by the organic species primarily results 
in the production of excitons rather than free electron-hole pairs 
{29]. These cells are cheaper to produce, lighter and more flexible 
than traditional silicon-based solar cells, and therefore, have a wider 
range of applications. 

Inorganic cells: In thin-film technologies, there exists a common 
problem with conversion efficiency due to poor material quality; 
photo-generated electrons and holes cannot travel very far before 
recombination (short free-carrier diffusion lengths). If a solar cell 
can be made using nanoscale hetero-junctions, then each photo- 
generated carrier will have less distance to travel, and the problem 
of recombination can be greatly reduced. The main factor for the 
superior efficiency of inorganic devices over organic devices lies in 
the high intrinsic carrier mobilities that exist in inorganic semi- 
conductors. Higher carrier mobility means that charges are trans- 
ported to the electrodes more quickly, which reduces current 
losses via recombination [30]. 

Quantum dot solar cells (QDSCs) are based on the Gratzel cell, or 
dye-sensitised solar cell, architecture but employ low band gap 
semiconductor nanoparticles that possess such small crystallite 
sizes that they form quantum dots (such as CdS, CdSe, Sb2S3 and 
PbS), instead of organic or organometallic dyes, as light absorbers. 
A quantum dot is a portion of matter that has electronic properties 
that are intermediate between those of bulk semiconductors and 
those of discrete molecules [31]. 


In general, silicon wafer-based photovoltaic cells that include 
single crystal, multi-crystalline and the hetero-structure have 
more than 80% of the market share [32]. The record conversion 
efficiency of a multi-crystalline photovoltaic cell is 20.4%, while a 
single crystal silicon cell can reach an efficiency of 26.4% (Table 1). 

To evaluate the performance and price of different PV module 
technologies above summarised, we can refer to the studies 
conducted by the International Energy Agency (IEA) exemplified 
in Fig. 1 [33], or in more recent studies [27]. 

The record conversion efficiency of amorphous silicon cells is 
13.4%, while the record efficiencies of CdTe and CulnSe> are 18.7% 
and 20.4%, respectively. Microcrystalline silicon has experienced a 
breakthrough recently and has reached an efficiency close to that 
of CdTe. Third-generation photovoltaic cells are solar cells that 
are potentially able to overcome the Shockley-Queisser limit of 
31-41% power efficiency for single band-gap solar cells. Common 
third-generation systems include multi-layer cells made of amor- 
phous silicon or gallium arsenide, whereas more theoretical 
developments include frequency conversion, hot-carrier effects 
and other multiple-carrier ejection. Third-generation photovoltaic 
cell technology is still in the early development stage. The cost- 
benefit ratio of this technology is still too high to be competitive 
with the previous generation counterpart. 


3. Generic electric behaviour of a solid single junction PV cell 


When a solid single junction PV cell is not illuminated, its 
behaviour is similar to a semiconductor junction; a simple diode, 
which has a characteristic I-V curve, is traditionally described by 
the equation: 


Ip = Io(ent? -1) (1) 


where I is the reverse saturation current of the diode, gis the 
electric charge of an electron (1.602 x 107"? C), kis the Boltzmann 
constant (1.381 x 107” J/K), T is the junction temperature, Vp is 
the voltage across the diode and n is the ideality factor, also known 
as the quality factor or sometimes the emission coefficient. The 
ideality factor n typically varies from 1 to 2 (though it can be 
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higher in some cases), depending on the fabrication process and 
semiconductor material; in many cases, n is assumed to be 
approximately equal to 1 (thus, the notation n often is omitted). 

The illumination of the semiconductor junction determines the 
rise of a photo-generated current J; that determines a vertically 
translation of the I-V curve of a quantity linked, almost exclusively, 
to the surface density of the incident energy (Fig. 2). 

In this way, an ideal cell is represented as a current generator, 
connected in parallel with a diode, and its I-V characteristic is 
described by Shockley [35] with the following equation: 


687 


Eq. (2) represents only a mere theoretical description because it 
does not take into account the effect due to the presence of the 
electrodes both above and below the semiconductor layer that are 
necessary for collecting the charges that partially cover the intercept- 
ing surface. In [36], Wolf correctly observed that in a PV cell, the 
photocurrent is generated not only by a single diode but also by a 
global effect due to the presence of a plurality of adjacent and 
uniformly distributed elementary diodes along the surface that 
separates the two layers of the semiconductor. Each elementary diode 
is traversed by a current that flows across the semiconductor layers 
along a different line, characterised by different electric resistance and 


Wp 
I=, -Ip =l —lo (ent 1) (2) voltage drop. 
here lis th it z The transverse component of the current l;, parallel to the 
WHhETe Toe Out UL CuETent. surface of the cell, should be different for each elementary diode; 
in this way, there is a different characteristic I-V for each 
Table 1 oe elementary diode. Because these diodes are considered in parallel, 
Best reseañch-cell efaendies of 2013/03184] the combination of them determines the global I-V characteristic 
Category Types Efficiency [%] of the PV cell. The transverse electrical resistance is much higher 
than the electrical resistance relative to the direct component of I;. 
Multi-junction cells Four junction or more 37.80 The transverse component of the current I,, which only occurs in a 
(2-teminal monolithic). (no concentrator) real photovoltaic cell, gives rise to a high energy dissipation that 
Three junction (concentrator) 44.00 d ically d h à ffici f th l Il 
Three junction (no concentrator) 37.70 rastica y decreases the conversion efficiency o t e solar cell. 
Two junction (concentrator) 32.60 Wolf, in [36], proposed a simplified equivalent circuit, shown in 
Two junction (no concentrator) 30.80 Fig. 3, in which there are only a couple of diodes, one current generator 
Single junction GaAs Concentrator 29.10 and two resistors that take account of the dissipative effects above 
Thin film crystal 28.80 described and of the presence of any constructive defects. 
Ser Ta E The resolution of the above equivalent circuit permits to obtain 
Single crysta 710 : . : x res i 
Thin film crystal 3010 the following expression of the I current in an implicit form: 
Crystalline Si cells Silicon heterostructures (HIT) 24.70 I= -I oar _1)-1 Ge rs 1) = V+IRs (3) 
Multi-crystalline 20.40 en 02 Rsh 
Thick Si film 16.00 à o 
Cu (In, Ge) Sez 20.40 where I; is proportional to the irradiance [35]; Io; and Io) are the 
Thin film technologies CdTe 18.70 diode saturation current of the two diodes, R, and R,,are the series 
Nano-micro-poly-Si 17.20 and parallel resistances, respectively; and n, and nz are the ideality 
Multi-junction poly-crystalline 15 factors of the two diodes. As it is easy to see from Fig. 3, it is necessary 
Amorphous SiH (stabilized) 13.40 to know seven parameters (1;,11, 112,101, lo2, Rs, Rsp) to solve Eq. (3). 
Emerging photovoltaic Dye-sensitized cells 11.40 In detail, the resistors R, and R, alter the slope of the curve 
emer 11.10 before and after the “knee” of the curve, respectively, while the 
as ziik D ratio between Ip, and loz changes the curvature [37]. 
Quantum dot cells 7.00 Despite not being a mathematically indeterminate problem, 
calculation of the seven parameters is made difficult by the 
USD USD USD USD USD Module 
200 /m? 300 /m? 400 /m? 500 /m? 600 /m? price per m? 
4L 
F 
3 3f 
= Crystalline 
a Silicon 
=} 5 
4 Technologies 
S 85-90% * 
E PEE Organic Thin Films 
Ej Solar Cells Technologies 
£ Under 1%* 10-15% * 
1 F 
(0) r - 7 7 > 
5% 10% 15% 20% 25% Efficiency 
50 W/m? 100 W/m? 150 W/m? 200 W/m? 250 W/m? Performance 


Fig. 1. Current performance and price of different PV module technologies [33]. 
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Fig. 3. Schema of simplified two diodes equivalent circuit. 


implicit form of the equation and by the presence of two 
exponential terms. In the scientific literature, there are few fully 
explained models that permit coding of the algorithm to retrieve 
the seven parameters; furthermore, these models concern the 
single cell and are based on particular assumptions that limit their 
application. Indeed, some authors admit that the initial conditions 
strongly affect the resolution of Eq. (3) [38-40]. In these models, 
the procedures need to be properly guided during the initial 
estimation of the parameters to avoid inconsistent results. 

For these reasons, some authors have preferred to employ 
different correlations not based on an electric model [41-42], 
whereas others have used a simplified model based on a single 
diode [15-19], as shown in Fig. 4. This new equivalent circuit is 
described by the following implicit equation: 


fei lofe" 1) (4) 


in which there are only five unknown parameters and only one 
exponential term. According to the traditional approach, the 
photocurrent I; depends on the irradiance, Ig is affected by the 
temperature of the cell and n, R; and R,, are constant. 

Some authors have focused on studying the one-diode model 
and have suggested some directions for improvement and/or 
simplification that would allow determination of the five para- 
meters (lz, 1,Jo, Rs, Rsp) based on the performance data of the 
modules commonly provided by the manufacturers. Here, we 
review five of the best explained models available in the literature 
for the evaluation of the electrical behaviour of a single junction- 
based PV module. The models proposed by De Blas et al. [15] 
(par. 3.1), Hadj Arab et al. [16] (par. 3.2), De Soto et al. [17] (par. 3.3), 
Villalva et al. [18] (par. 3.4) and Lo Brano et al. | 19] (par. 3.5) will be 
described in detail. 


Fig. 4. Schema of one-diode simplified equivalent circuit. 


3.1. De Blas, Torres, Pieto and Garcia model 


In [15], De Blas et al. address the study of the electrical 
behaviour of a PV panel determining the five parameters solving 
an equations system and applying some modifications and sim- 
plifications to the traditional model. De Blas et al. state that their 
model permits to obtain a series of I-V curves taking into account 
the environmental variability with errors less than 1% in compar- 
ison with experimental results. 

Three of the five equations of the system, are determined 
imposing that the following points belong to the I-V curve, 
expressed by Eq. (4): 


1. Open circuit point: (V = Voc, I= 0); 
2. Short circuit point (V = 0, I= Isc); 


3. Maximum power point(V = Vmp, [= Imp). 


In these cases, Eq. (4) becomes 


Omi lo (eF 1) = (5) 
S 
Ise =I (e4* 1) -55 6) 
S| 
| eens a (ene ve 1) - fuels mpRs (7) 
S| 


The substitution in Eq. (6) of the value of the photocurrent derived 
from Eq. (5) yields the following: 


Ise(14+g¢) <1 lo (e° 1) =1o(e* ja lo (e 1) 


Rsn 
Rs Voc Voc IscRs 
Ise(1+3") =h(e E1 e +1) 


Then, 


„Rs \ Voc] ,/ ¥oc ses 
e- bee) Ha e-e 


assuming that 


V, TecR 
et sen (8) 


De Blas et al. hypothesise that 


_ Rs Voc] , “oc 

” E isc(1 A Rsh ie i © 
Similarly, substituting Eq. (7) the same value of the photocurrent, 
the equation becomes 


Yoc _Yoc Vmp + tmp Rs Vinp + ImpRs 
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Furthermore, substituting Eq. (9) in Eq. (10) yields the following: 


Rs oa Vp + ImpRs ~ Voc Voc— Vinp ImpRs 
Imp = |I 14g) - {er |g OC mp ns 
me | sE ( Rsh Rsh | | Rsh Rsn 


or else 


Rs ) ( Rs ) Voc Voc—Vmp 
Imp| 1+ Isc{ 1+ + 
ie ( Rsh j Rsn Rsh Rsn 


Rs Voc Vp +ImpRs -Voc 
E Isc(1+R2) Rsh |! E j | 


After some mathematical manipulation, it is possible to state the 
following: 


Vinp +InpRs = Voc 


= [np (Cre Z590 + Rs /Rsn) = (V mp /Rsn) 
Isc(1 +Rs/Rsn) — (Voc /Rsn) 


nT 


(11) 


The other two conditions that permit solving the equation system 
are determined by imposing the equality of the derivative of 
Eq. (4) at short circuit and open circuit points, respectively, with 
the reciprocal of the slopes of the characteristic curve in these 
points; these values are called Rmo and Ro, and they can be 
graphically evaluated from the experimental I-V curves issued by 
the manufacturers. 


4. The derivative at the short circuit point is as follows: 


ls 


A ps 1 
dI nteXP "™ +R 1 


dV 


(12) 


E = Io IscRs 1\ R 
V=0 1+R, (exp nT +R shO 


5. The derivative at the open circuit point is as follows: 


dl mete 1 aa 
AV y = voc 1+R, (Hem +) Rs, 


3.1.1. Mathematical and physical assumptions 

Furthermore, to simplify the equation system, De Blas hypothe- 
sises that 
Io IscRs 1 


ew < 
nT sh 


(14) 


so Eq. (12) becomes 

dI 1/R;n 1 

dV|lyvz0 1+(Rs/Rsn) Rho’ 

Finally, it is possible to state that 

Rsn = Rsho — Rs (15) 


Substituting Eq. (9) into Eq. (13) yields the following: 


Isc( Ron + Rs) ~Yoc 


and because the value of Eq. (15) can be expressed as function of 
the resistances Rmo and Ro, the last equation can be rewritten as 
follows: 


1 Isc(Rsho =R HR) -Voc +1 


Rs, Rsno — Rs + Rs [Rao Roe Rs) —Voc 4. 1] 


IscRsno — V IscRsno — Vi 
R£ Ha OC E Runo = Ro [r+ oa =. 


Explicating R, results in the following: 


p Pal tat + Rt + 


Voc = IscRsno 
nT 


(16) 


The previous simplifications that have been adopted to determine 
the unknown parameters, e**>e*r and ‘er <; are usually 
sh 


n 
fully satisfied. 


3.1.2. The algorithm steps 

With this approach, adopting an iterative process that provides 
the initial estimate of R, value, the De Blas et al. model retrieves 
the five unknown parameters. In detail, the algorithm steps are as 
follows: 


a starting R, value is hypothesised; 

n and Rs, values are calculated from Eqs. (11) and (15); 

Rs value is recalculated with Eq. (16); 

an iterative procedure verifies the convergence of R, value; 
when R, converges, the correct value of R,, is recalculated; and 
the values of I, and Ip.are calculated using Eqs. (5) and (9). 


To describe the behaviour of a PV module at different condi- 
tions of temperature and radiation, De Blas et al. suggest the 
application of the procedure described by the International Stan- 
dards IEC891 [43] that, starting from the knowledge of current and 
voltage values for the given temperature and irradiance, defines 
the variation of the I-V characteristics. To adapt the saturation 
current of the diode Jo with a variable temperature conditions, 
IEC891 uses the Townsend equation [44]: 


I(T) =D Tet) (17) 


in which €;e is the material band-gap energy (1.12 eV for silicon ) 
at Standard Reference Conditions (SRC), A=y N, where y is the 
completion or shape factor and N, is the number of cells wired in 
series (for ideal cell A is equal 1) and D is the diode diffusion factor 
(approximately constant). Therefore, on the basis of the correct 
evaluation oflọ, Rs and R,,, the other parameters can be derived 
using the values of Voc and Isc at new conditions (Eqs. (5) and (6)). 

With these assumptions, the parameters that are variable with 
temperature and/or irradiance are lo, Isc and n. This last parameter, 
which should be considered independent of temperature, can be 
evaluated only if the values of Imp and Vmp are known in non- 
standard conditions. 


3.2. Hadj Arab, Chenlo and Benghanem model 


In [16], Hadj Arab et al. proposed a five parameter model that is 
detailed described and experimented in [45]. In this last work the 
electrical behaviour of a PV panel by referring to an equivalent 
electric circuit which includes the photocurrent, the reverse 
saturation current and the series and in parallel resistances. 
Indeed the experimental test show how the traditional model 
with five parameters leads to more accurate assessments of 
operating current and more faithfully describes the behaviour of 
a photovoltaic panel especially during the maximum irradiance 

In this case, Eq. (4) has been rewritten in a different form: 


V+IR; 
Rsh 


faa, (Gag 1) (18) 


in which the five unknown parameters of the model proposed in 
[25] are I;, Io, Rs, Rn and m when 
[Vinp +ImpRso = Voc] 
g Vm V Im 
Ve [In(Isc Rn Imp) In (Isc re) +w] 


a9) 
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To determine the values of these parameters, the three known I-V 
pairs at SRC are imposed in Eq. (18): open circuit point 
(V= Voc, I= 0); short circuit point (V=0, I = Isc) and maximum 
power point (V = Vmp, I= Imp). Furthermore the resistance defini- 
tions used in Eqs. (12) and (13) are applied; where Ro and Rṣno are 
the inverse of the slope in the open circuit point and in the short- 
circuit point, respectively; these assumptions lead to a different 
form of Eqs. (12) and (13) in which nT is substituted by mV, 


3.2.1. Mathematical and physical assumptions 
To simplify the solution of this system in the model the 
following considerations are applied: 


1 gga (20) 


Eq. (13) becomes 


b ptt 
m 
mE 5 1 


b oe 
1+ Ropper 


e: =0 22 


Eq. (12) can be rewritten as follows: 


1/Rsn 1. 1 1 
1+ (Rs/Rsh) z Rsho’ Rs +Rsh E Rsno 


So, 
Rsh = Rsho —Rs (23) 


Furthermore if 

R¢<Rsh > Rsho & Rsh (24) 
Similarly, if 

er>] 25) 
in Eq. (5), and I; = Isc, then 


I, —Yoc 


V 
Ip =—* = (Isc neem (26) 


Knowing Io, Eq. (7) can be rewritten as follows: 


v, 
Isc — 2° (Vp + Imps v LR 
Imp =lÍsc— Wee e ™ —1)/—-=" re z 


emv; 
V Vp + ImpRs — Voc Voc 

=lsc— (Isc OC) (e om eTe 
Rsh Rsh 


Vmp + ImpRs 


e m >] 27 


If 


V, 
then the term e~™ can be neglected: 


Voc Vp + ImpRs — Voc Vinp Rs 
Isc — me =! I 1 
( z Rsh ) ° = Rsh ki Rsh 


Finally, if 


Rs 
Rsh 


<1 (28) 


using the definition of Rs in Eq. (21), 


ImpmvV, 
Vata clip sy = TE 


v ce v 
lege | ae E ngs 
SC R =45C R mp> 
sh sh 
imp v 
Vmp +ImpRso -Voc ç Yor (Isc— 32- Imp) 
e mVt e Rsh = he — 
Oc 
(Isc = vec) 
or rather 
Isc- ee —I 
Ving timp R-Voc -= ImpmV: ia | sc- Re =| (29) 
V, ie — Voc 
mVt mV; (Is a) Isc Ren 


Thus, it is possible to obtain the definition of mV; as explained in 
Eq. (19). 

To solve the system formed by Eqs. (6), (19), (24), (29) and (26), 
is necessary to know the values of Isc and Voc. These parameters 
are considered variables with irradiance and temperature: 


G 
Isc = Ise.src Cac + M),.(T —T src) (30) 
G 
Voc = Voc.src + mV; In Gear + Hvac (T = Tsrc) (31) 
SRC 


3.2.2. The algorithm steps 
The calculation procedure can be summarised as follows: 


© the R,, value is calculated by Eq. (24); 

Isc and Voc parameters are calculated by Eqs. (30) and (31); 

e@ mV; is obtained knowing Ryo , Isc, Voc, Current and voltage at 
the maximum power point and using Eq. (19); 

© R, is obtained with Eq. (21); and 

Io is obtained from Eq. (26); 

è Finally, I, is obtained from Eq. (6). 


3.3. De Soto, Klein and Beckman model 


De Soto et al. propose a model in [17] for the evaluation of the 
performance of a PV panel represented by a one-diode equivalent 
circuit in which the current-voltage relationship for fixed cell 
temperature and irradiance is expressed by the following equation: 


VIR V+IRs 
=h (e z 1) (32) 
where a is an ideality factor defined as follows: 
iz Se (33) 


Even in this model, the solution of the five parameter problem 
implies the application of some conditions that refer to the points 
of open circuit (V = Voc, I= 0), short circuit (V = 0, I= Isc)and the 
current and voltage at the maximum power point (V = Vmp, 
I= Imp) evaluated at SRC. The first three equations of the system 
are as follows: 


I= (e 1) K =0 (34) 
S| 
kesi =In(e*r -3 1) z “oe (35) 
S| 
B Vmp + Imp Rs Vinp +Imp Rs 
Imp = I, —lo(€ -1) = Ra (36) 


where a represents the ideality factor and I; is the value of the 
photocurrent under SRC. 
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Because these relations are insufficient for the determination 
of the five parameters (l, n, Io, Rs, Rsn), further equations are 
obtained by imposing the derivative of the power curve at the 
maximum point equal to zero and giving as known the coefficient 
value of the thermal variation of Voc: 


dP dI 
adv PaPa =Imp + Vmpay 


=0 (37) 
mp 


The last conditions use the expression of the thermal drift of the 
open-circuit voltage: 


dV 


V -V 
Hva = FF ~ Y OC, SRC — V OCT (38) 


1=0 Tsrc—T 


where Tsrc is the temperature at SRC. 

To determine jy, it is necessary to know the value of the open 
circuit voltage Voc. Because the parallel resistance R,, was 
assumed to be independent of temperature, to determine the 
value of Vocr using Eq. (32) and then to solve Eq. (37), it is 
necessary to know the dependency of the parameters a, Ip and Iz, 
with temperature. 


3.3.1. Mathematical and physical assumptions 
To this aim, the De Soto model takes into account the following 
assumptions: 


a) The value of ideality factor a is deduced from the temperature 
of the cell T. In particular, for any operating condition, it is a 
linear function of the temperature of the cell: 

a T 
== 39 
asrc Tsrc G3) 


b) The variation of the current with the operating temperature is 
explained by the following equation: 


akee 
Ío <| T als- n 
Tosrc  (Tsrc 


in which the ratio between the current at any conditions and 
the current obtained at reference conditions Ig sęc depends on 
the ratio between the temperatures and on an exponential 
term. 

In detail, the reverse current saturation is evaluated with a 
relation that is derived from the Townsend equation Eq. (17). 
In Millman-Halkias [46], it is stated that 


Ip = K Ten (41) 


where K is a constant and qVco = €src is the amplitude in joule 
of energy gap; for silicon, the values of 7 and m are 2 and 1.5, 
respectively. Substituting the value of Vr = kTN,/q into Eq. (17) 
and Eq. (41), it is possible to state that K =D, m=3,y=n. 
Furthermore, for an ideal cell, the equation becomes 


3 ope 
Tosrc =D Tspce* Tsrc (42) 
Substituting Eq. (42) into Eq. (40), it is possible to obtain the 


relationship that describes the dependence of the current Jp on 
the operating temperature of the system: 


ek | Tre 


378 3 alee 7 
b DT2ett | a EE (43) 


= — 
losrce prera — LT src 


c) The determination of the photocurrent I; is obtained as a 
function of short-circuit current thermal drift y;,.. In particular, 


the expression that allows the determination of this variable is 
as follows: 


Ir src +H T- Tsro)] (44) 


I -< 
J Gsrc Msrc 


where G represents the irradiance and M is the modified air mass 
(M is an index derived as function of air mass, zenith angle of the 
photovoltaic system and of four constants typical of each type of 
panel; these parameters were experimentally obtained from NIST 
- Fanney et al. [47]). 

Finally, to complete the model, De Soto et al. have analysed the 
dependence of the resistances R, and R,, with the operating 
conditions. The R, influences the shape of the I-V curves near the 
maximum power point. In [17], De Soto et al. show how, by 
varying R, by 20% (in excess or in deficit) with respect to the 
reference value, the effect on the curve is practically insignificant, 
thus justifying the assumption that R, can be treated as a constant. 
R;, determines the slope of the I-V curve in the short circuit 
condition, becoming more horizontal as it increases. In [17], the 
authors claim that its value seems to vary with the irradiance G for 
all types of cells; in particular, De Soto et al. refer to the Schroder 
definition [48], which indicates that, for low values of irradiance, 
the value of R,, is inversely proportional to the short-circuit 
current and therefore to the irradiance. For a description of this 
phenomenon, the De Soto et al. model proposes the following 
expression: 


Rsn Gsrc 
= ESRC 45 
Rsnsrc G (9) 


3.3.2. The algorithm steps 

Knowing the dependence of a, Ig and J, on temperature 
(Eqs. (39), (40) and (44)), and of R; and R,, on the operating 
conditions, the solution of the system of five equations (Eqs. (34)-(37) 
and (38)) with five unknowns is solved by a nonlinear solver, the 
Engineering Equation Solver — EES [49]. 

However, the operating procedure adopted to solve this equa- 
tions system represents a weak point in the model. 


3.4. Villalva, Gazoli and Filho model 


In [18], Villalva et al. propose a model to determine the 
parameters of the nonlinear I-V, shown in Eq. (32), by adjusting 
the curve at three remarkable points: open circuit, maximum 
power and short circuit points provided by the manufacturer. 
Given these three points, which are provided in all commercial 
array datasheets, the method aims to find the best J-V equation 
including the effect of the series and parallel resistances, and 
warranties that the maximum power of the model matches the 
maximum power of the real array. 


3.4.1. Mathematical and physical assumptions 
Two other conditions are necessary to solve the system of 
equations. In [18,50], the authors propose the following assumptions: 


© The I, current of the PV cell depends linearly on the irradiance 
and the temperature according to the following equation: 


I = [IL src + Hi (T— Tsg (46) 
© The authors propose a new expression of Io : 
I T-T 
Io = SEC +H, SRC) (47) 


v + T-T, 
OC,SRC +V oc (T ~Tsrc) 
e Ve -1 
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Usually, 1 <a < 1.5 because a expresses the degree of ideality 
of the diode, and it is totally empirical; any initial value can be 
chosen to adjust the model. Eq. (47) proposes a different approach 
to express the dependence of Jo on temperature and the thermal 
drift of voltage. 


3.4.2. The algorithm steps 

In [18,50], the authors propose an algorithm for adjusting the 
resistance values based on the fact that there is only one pair {R,, 
Rsn} that assures that the calculated maximum power is equal to 
the maximum experimental power. The authors claim that the 
proposed model is further improved by taking advantage of the 
following iterative solution of R, and Rsn: Rs is slowly incremented 
starting from R;~0; the initial value of R,, is given by the 
following: 


Vinp Vocsrc — Vinp (48) 
Iscsrc = Imp Imp 


Rsh, min = 
Each iteration updates R, and R,, towards the best model solution. 
The initial guess of I; src is obtained by the following relationship: 


I, pee Rs + Rsn 
Rsh 


Ísc,src (49) 


The values of R, and R,, are initially unknown, but as the solution 
of the algorithm is refined along successive iterations, the values of 
R, and R,n tend towards the best solution, and Eq. (49) becomes 
valid and effectively determines the light-generated current I; 
while taking into account the influence of the series and parallel 
resistances of the array. The software coding this algorithm is 
available in [51]. 


3.5. Lo Brano, Orioli, Ciulla and Di Gangi model 


Lo Brano et al. propose, in [19], a model for the evaluation of 
the performance of a PV panel represented by a one-diode 
equivalent circuit in which it is postulated that Rs, R,, and Io 
change with the irradiance; this hypothesis is realistic because a 
PV cell and an illuminated semiconductor diode have different 
electrical behaviours. Based on these assumptions, a new five 
parameter model is proposed, and it is described by the following 
equation: 


ac[V +J I(T-Tsrc)] +! Rs 
Kac, T) = acli(T)— lo(ac, T) (e agi = ) 
ac[V +] I(T —Tsrc)| +1 Rs 


Rsh 60 


where a@cg=G/Gsrc represents the ratio between the generic 
irradiance G and the irradiance at SRC Gsrc, and I(T) is the 
photocurrent that can be evaluated by the following relation: 


I(T) = IL sre + My,.(T — Tsro) (51) 


3.5.1. Mathematical and physical assumptions 

Although the authors suggest a linear relationship between the 
values of the series and shunt resistances and the irradiance, the 
following positions are made: 


ne 


52 
Ron (ac) = “ast se 


Another important innovation concerns the determination of 
the value of Io(ac¢), which can be represented by the following 


equation: 


i= Voc(&c,T) 
iac, T= a6“ (53) 
em 
The variability of Voc with the temperature can be evaluated with 


the following equation: 
Voc(a&c, T) = Voc,src(&c, T) + My,.(T —T src) (54) 


Lo Brano et al. observed, in [19], that for numerous commercial 
panels, there is a regular dependence of the reverse saturation 
current with the irradiance and temperature expressed by the 
following formula: 


ac — 200) 
1202 


Ig(1:T) 
™79(0.2:T) 


Io(æg; T) = exp ( + In Io(0.2; | (55) 
To better reproduce the I-V characteristic at generic temperature T, a 
thermal correction factor J is defined as follows: 


Vinpsrc — Vmp(&G, T) 
a Imp(T — T src) ee 


Here J is a thermal correction factor similar to the fill factor of the 
curve described by the IEC891 [43]. To determine the value of J, the 
value of T should be chosen based on the maximum or the 
minimum working temperature of the PV module and the max- 
imum range of temperature within which the data from the 
manufacturer are provided (in [19], a temperature of 75 °C is 
considered). 


3.5.2. The algorithm steps 

The procedure to determine the five parameters requires a 
preliminary assessment of R; and R,;,. It is also necessary to known 
the values of Isc.srcs Voc.src lmp,sRC> V mp,SRC and the value of Voc at 
standard temperature and at the lowest irradiance (typically 
200 W/m°); all these values can be graphically extracted from 
the standard I-V characteristic curve provided by the manufac- 
turer. To determine J it is necessary to have the values of Imp and 
Vmp at different temperatures with respect to SRC. To determine 
the values of the unknown parameters, the three known I-V pairs 
at SRC are imposed on Eq. (50), obtaining the same functional 
relationships as the Blas et al. model (Eqs. (5)-(7), (12) and (13)). 

The system is solved numerically via a double cycle of iterations 
(trial and error) to retrieve the values of R, and n. To initiate the 
solving procedure, the following positions are made: 


l Ir src = Isc.src 


57 
Rsh = Rsho 6% 


which is obtained by simplifying Eqs. (13) and (5) and taking into 
account the following: 
Rs<Rsho 
TosrcRs (58) 


T. 
Tosrce "SRC 1 
nT sxc Rsno 


The procedure comprises the following steps: 


1. first attempt values of n and of R, are set; 

. ILsrc and R,, are calculated with Eq. (57); 

3. Tosrc, Insrc and R,, are derived in sequence, respectively from 
Eqs. (7), (6), and (12); 

4. n is derived from Eq. (5), compared with the previous attempt 
value and is changed accordingly; 

5. the trial and error process in steps 3 and 4 is reiterated until 
convergence of n is achieved with the desired accuracy; 

6. once the first trial and error process for n has been completed, 
Rs is derived from Eq. (12), compared with the previous 
attempt value and changed accordingly; 


N 
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Table 2 
Synopsis of mathematical relationship employed for each model. 


Five parameter models 


De Blas et al. 


I-V relationship for 1st, 2nd and 3rd conditions 


4th condition 5th condition 


; Vte V+IRs iri ig Yop 
Hadj Arab et al. T=I,—Io (e r 1) Rar m _ exp T rh al E Ben +r 
P _ dviv=o0— 7 ( Tc Rs ) WIV =Vo ( Yoc ) 
T=acl,—Io eat 1) 1+R; exp nT ty of 1+Rs Denr + 
1 1 
Lo Brano et al. I=I(sc,0) =0 ac[V + JI(T —Tspc)]+IRs as =- 
si 0 
1=1(0,Voc) =0 Rsh 
1=I(Imp, Vp) =0 
dP = — dv a Vocsrc —Vocr 
De Soto et al. leh p (e1) YeR IP = Prox = Moe = Gli =0 ~ SE 
i Emoe ) Kn Isesre + Hig (T — Tse) 
Villalva et al. ( sn i=j T-T G sc src + Pige (T— Tec 
L = [lL src + Mice ( RONG lowa tiv Tae Vy cao Fan Teac WV A 


7. the first trial and error process in steps 3-5 is again started and 
completely developed; 

8. Rs is again derived from Eq. (12) and compared with the value 
of the previous attempt and changed accordingly; and 

9. the second trial and error process for R; is reiterated until the 
desired accuracy is achieved. 


The system is solved in a non-approximate form because no 
simplified equations were used. The above procedure is coded into 
a short software routine using simple languages. It does not need 
to be guided towards solutions starting from initial fitted values of 
the searched parameters. The procedure is more reliable because it 
normalises all parameters with respect to Voc and Isc. Further- 
more, in [52], the Lo Brano et al. model was validated with an 
experimental setup. 


4. Discussion 


In this section, the models described above are compared, 
highlighting their peculiarities while discussing the availability 
of initial data and the reliability of the results. 


4.1. Availability of necessary input data 


Each model is characterised by a system of five equations with 
five unknowns (I, n, Io, Rs and Rsn). All considered models require 
the following three known points of the I-V curve: open circuit 
point, short circuit point and maximum power point. These points 
are always available in the datasheets provided by the manufac- 
turers, and therefore, the obtainment of these necessary values is 
immediate. 


è The De Blas, Hadj Arab and Lo Brano models require the 
measurement of the slope of the I-V curve at selected points: 
open and short circuit. Normally, the technical datasheets 
provided by the manufacturer do not include these data, and 
a graphic valuation must be carried out. This procedure is only 
possible if the datasheets provide a graphical representation of 
the I-V curve; however, a graphical estimation is always 
approximate, which produces a measurement error that can 
affect the results: 

© The De Soto model also needs a graphical estimation of the I-V 
slope at the maximum power point and also requires the value 
of Voc, the thermal drift coefficient. Although this value is 
generally available, relationship (38) is a finite numerical 
approximation of a derivative. 

© On the other hand, the Villalva model requires the thermal drift 
coefficients of the open circuit voltage and short circuit current. 
These values are generally available in the official datasheets of 


Table 3 
Mathematical and physical assumptions for each model. 


Assumptions 


De Blas et al. Eqs. (8), (14), (17) 


Hadj Arab et al. Eqs. (20), (22), (24), (25), (27), (28), (30), (31) 
Lo Brano et al. Eqs. (52)-(58) 

De Soto et al. Eqs. (39), (43), (44) 

Villalva et al. Eqs. (46)-(49) 


commercial PV modules. This model does not require any 
graphical estimation. 


4.2. Synopsis of the major hypothesis and simplifications 


To facilitate the identification of the mathematical relationships 
that permit coding of the models into a software package, 
Tables 2 and 3 collect, for each model, the five fundamental 
equations used and their simplifications. As summarised in the 
first column of Table 2, all the models employ three remarkable 
known l-V pairs at SRC. These values are substituted in the 
fundamental equation; the second column of Table 2 collects for 
each model the mathematical function that links current and 
voltage, I= I(I, V). The last two columns gather the fourth and 
the fifth conditions, which permits solution of the equation 
system: 


è The De Blas et al. and the Hadj Arab et al. models are 
characterised by the same system of equations; however, the 
different physical and mathematical assumptions listed in 
Table 3 lead to different results. The fourth and fifth conditions 
concern the graphical valuation of I-V slope at short circuit and 
open circuit points. 

© The first three equations employed by the De Soto at al. and 
Villalva et al. models are identical; nevertheless, the models use 
the thermal drift coefficients in different ways. The De Soto 
et al. model also uses the slope of the P-V curve at the 
maximum power point. 

© Finally, the Lo Brano et al. model proposes a reformulation of 
the general fundamental equation that directly takes into 
account the effect of temperature and irradiance. The last two 
conditions are identical to those used by the De Blas et al. and 
the Hadj Arab et al. models. 


4.3. Numerical performance 
To assess the accuracy of each model with respect to the data 


provided by the manufacturer, the I-V and P-V curves at vari- 
able temperature and irradiance obtained by each model are 
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Fig. 6. I-V curves of a PV panel at variable irradiance and constant temperature. 
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Fig. 7. P-V curves of a PV panel at variable irradiance and constant temperature. 


represented in the following graphs (Figs. 6-9). In more detail, the 
authors used a specimen PV panel with the I-V curves drawn in 
Fig. 5 and the tabular data collected in Table 4. 

Due to a lack of information regarding a detailed calculation 
procedure, it was impossible to obtain any numerical results from 
the De Soto et al. model. 

The graphs above show that all of the models describe the I-V 
and P-V curves well when the irradiance is 1000 W/m? and the 
cell temperature is 25 °C. By varying the irradiance and/or the 
temperature, the agreement with the curves provided by the 
manufacturer changes considerably. Defining the mean absolute 
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Fig. 8. I-V curves of PV panel at variable temperature and constant irradiance. 
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Fig. 9. P-V curves of a PV panel at variable temperature and constant irradiance. 


error as follows (Tables 5 and 6): 


ye P; calculated —P; manufacturer 
MAE = : : 
N 


i=1 P; manufacturer 


where N is the number of the samples, the following tables show 
the accuracy of each model with respect to the data issued by the 
manufacturer: 


© Generally, the model that best describes the I-V and the P-V 
curves provided by the manufacturer is the Lo Brano et al. 
model, which is characterised by an MAE of 0.59 for the power 


G. Ciulla et al. / Renewable and Sustainable Energy Reviews 32 (2014) 684-696 695 


Table 4 
Specimen PV panel characteristics. 


Voc [V] Isc [A] Vmax [V] Imax [A] #voc [V/°C] Hise [A[°C] 


1000 W/m? 29.35 8.07 23.6 7.567 —1.07 x 107! 2.22 x 107? 
200 W/m? 27.2 1.614 22 1.508 


Table 5 
Mean absolute error at constant temperature (25 °C) between power data issued by 
the manufacturer and each model results. 


Mean average error Global MAE 
200 400 600 800 1000 
W/m? K Wm?K Wm?K Wm?K Wm? K 
De Blas 1.36 1.87 2.02 1.48 0.39 1.42 
Hadj Arab 3.19 4.82 5.01 4.12 0.39 3.51 
Lo Brano 0.41 0.68 1.06 0.40 0.39 0.59 
Villalva 1.18 2.33 2.38 2.21 1.70 1.96 


Table 6 
Mean absolute error at constant irradiance (1000 W/m?) between power data 
issued by the manufacturer and each model results. 


Mean average error Global MAE 
25°C 50°C 75°C 
Hadj Arab 0.39 2.07 3.71 2.05 
Lo Brano 0.39 0.53 0.82 0.58 
Villalva 1.70 1.16 1.96 1.61 


calculated at different values of irradiance and constant tem- 
perature and by an MAE of 0.58 for constant irradiance and 
different temperatures. 

© The De Blas et al. and Hadj Arab et al. models tend to slightly 
underestimate the current and power values before the max- 
imum power point; however, after the maximum power con- 
ditions, both models overestimate the currents and shift the 
curve slightly to the right. The De Blas model does not allow for 
calculation of the curve at variable temperatures. 

è At the same time, the Villalva et al. model tends to overestimate 
the current and power values under any conditions other 
than SRC. 


4.4. Overview of the models, hypothesis and simplifications 


All examined models employed an iterative procedure to solve 
a system of five equations. However, each model is characterised 
by specific assumptions and simplifications that inevitably affect 
the results: 


© The De Blas et al. model postulates that the quality factor n, 
which should be constant, varies with temperature and irra- 
diance; in general, the value of n is known only at SRC (Eq. 
(11)). In practice, for each condition of irradiance and tempera- 
ture, the model needs to be recalculated, estimating again the 
values of Ro and R,pg. To avoid uncertainty in the application of 
corrective methods, in [15], De Blas et al. validated their model 
using Isc and Voc values known under specific conditions of 
temperature and irradiance. 

è The Hadj Arab et al. model is tested in [16]. This latest work is 
the one that best describes the proposed model and that best 


matches the electrical behaviour of a solar panel. Nevertheless, 
the model is characterised by several simplifications regarding 
the exponential terms that could affect the results. Further- 
more, the model requires that the shunt resistance be much 
greater than the series resistance. This condition is generally 
always verified in a high quality PV module. 

© In De Soto et al., the five parameters are determined through a 
nonlinear solver, a calculation tool that may have difficulty in 
solving the equation system. Some tests performed with auto- 
matic solvers have not led to automatic convergence. In 
addition, in the original article [17], there is not enough 
information about the calculation steps. 

© Lo Brano et al. proposed a new five parameter model that was 
improved and corrected. The model is characterised by some 
significant differences, such as the variation of Ip, Rs and Rs, 
with irradiance and temperature. Application of the model 
allows for obtaining the best results in terms of comparison 
with the experimental I-V and P-V curves issued by 
manufactures. 

è The aim of the Villalva at al. model is to fit the mathematical I- 
V equation to the experimentally remarkable point of the I-V 
curve of the practical array. A new empirical equation is 
proposed to calculate the reverse diode saturation current. 
Furthermore, an algorithm coded into a routine is available 
on the web for public download. 


5. Concluding remarks 


The I-V characteristics of a photovoltaic PV module can be 
reproduced assuming a one-diode equivalent circuit made of 
linear and non-linear components. In this paper, five of the most 
recent mathematical one-diode models for predicting the I-V 
curve of a PV cell, which take into account the variation of 
irradiance, electrical load and operating temperature, were exam- 
ined in detail. These models reduce the problem to the determina- 
tion of only five parameters: the photocurrent I, the diode 
saturation current Ip, the ideality factor n, and the series and 
parallel resistances R, and Rgp. 

The procedure for calculating these five electrical parameters 
constitutes the heart of the models examined, and this work 
allows for their evaluation with regard to the fundamental 
hypothesis, the availability of necessary data and the reliability 
of the results. For each model, the fundamental equations and the 
most important assumptions used were reported, and to simplify 
the analysis of the procedures, all equations were rewritten with 
the same nomenclature where possible. Finally, an accuracy 
analysis was performed to assess the capability of each model to 
reproduce the I-V curve of a specimen PV panel. This work 
permitted verification that all the models considered require 
numerical data that are not always available in the official 
datasheets. Some models require data obtainable only by applying 
semi-graphic procedures, thus leading to an inevitable uncertainty 
in the measured values. A common feature of all models evaluated 
in this paper is that, for solution of the system of five equations, 
they require the application of numerical techniques that may not 
converge if the first guess values are not sufficiently close to the 
real solution. 

Concerning the reliability of the results, all of the models can 
reproduce the I-V curve at SRC with sufficient accuracy, although 
the Lo Brano et al. model is slightly more accurate than the others. 
However, although it is not particularly accurate in the calculation 
of the I-V and P-V curves, the Villalva et al. model does not require 
any graphic measurements, and a script is already available on the 
web. 
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